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Abstract 

Using the Lorentz Integral Transform (LIT) method we compare the results for the triton total 
photodisintegration cross section obtained using the Correlated Hyper spherical Harmonics (CHH) 
and the Effective Interaction Hyper spherical Harmonics (EIHH) techniques. We show that these 
two approaches, while rather different both conceptually and computationally, lead to results which 
coincide within high accuracy. The calculations which include two- and three-body forces are of 
the same high quality in both cases. We also discuss the comparison of the two approaches in 
terms of computational efficiency. These results are of major importance in view of applications 
to the much debated case of the four-nucleon photoabsorption. 

PACS numbers: 21.45.+V, 25.10.+S, 25.20.DC, 27.10.+h 
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I. INTRODUCTION 



In the last few years there has be a resurgence of interest in the photodisintegration of 
few-nucleon systems. This is essentially due to two reasons. One of them is that few-body 
theories have now achieved the ability to produce results of high accuracy for reactions 
into the many-body continuum. The other is the ongoing debate about the existence and 
nature of three-body forces, which finds a natural testing ground in reactions involving the 
three-nucleon systems. Some theoretical results on photonuclear reactions with light nuclei 
using realistic two- and three-body potentials have already been published [? ? ? ] and 
among them is a benchmark between the correlated hyperspherical harmonics (CHH) and 
Faddeev approaches [? ] on the total photodisintegration cross section of the three-body 
systems. This benchmark demonstrated the high degree of accuracy reached by these two 
techniques. At the same time it showed the inadequacy of existing experimental data to 
allow any conclusion about the three-nucleon force effect, whose magnitude turns out to be 
smaller than the rather large error bars. 

In this work we test the EIHH method on the triton photonuclear cross section with 
realistic forces. This method is also based on a hyperspherical harmonics expansion, but 
uses the concept of an effective interaction to speed up the convergence of the expansion. 
We believe that it is necessary to benchmark EIHH results on a realistic force calculation. 
In fact, when combined with the LIT method [? ], this approach is very promising (see e.g. 
[? ? ]) for the study of electromagnetic reactions on heavier systems, where three-body 
force effects might be more prominent. 

In the following we first briefly recall the LIT method (Section IHJ. In Section ITTTl the 
hyperspherical expansion method is summarized, together with the CHH fSection llll A|) and 
EIHH f Section IIII B|) techniques. They address the problem of convergence for the solutions 
of the Schrodinger and Schrodinger-like equations required by the LIT method. Section 
IIVI contains results for the triton photoabsorption with a realistic two- and three-nucleon 
force and, in view of applications to the electromagnetic responses of heavier systems, also 
a discussion about the computational efficiency of the two approaches. 



3 



II. THE LIT METHOD 



The LIT method has been described extensively in several publications [? ? ]• Here we 
only recall the key points. 

The calculation of any inclusive electromagnetic cross section requires the knowledge of 
the response function 

Ri^) = f^.m^M'Siu -E^ + Eo) , (1) 



where 00 represents the energy transferred by the elctromagnetic probe, |0) and Eq are ground 
state wave function and energy of the system undergoing the reaction, \n) and En denote 
eigenstates and eigenvalues of the Hamiltonian H, and G is the operator relevant to the 
reaction. The relation between the total photoabsorption cross section and the response 
function is given by 

= 47i^auR{uj) , (2) 

where a is the fine structure constant and G of Eq. (PJ) is the nuclear current operator. In 
the low-energy region considered here one can rely on Siegert's theorem and use for G the 
unretarded dipole operator 

A 

e = J2^^rf. (3) 

i 

The LIT method consists in calculating R{uj) in three steps. Step 1. The equation 

{H~E,-uj, + iT)\^>)=Q\Q). (4) 

is solved for many and a fixed F. 

This is a Schrodinger-like equation with a source. It can be shown easily that the solutions 
are localized. Thus one only needs a bound state technique to obtain them. In particular 
one can adopt the same bound state technique as for the solution of the ground state, which 
is an input for Eq. (jlj). 

The values of the parameters ujq and F are chosen in relation to the physical problem. 
In fact, as will become clear in Step 2, the value of F is a kind of "energy resolution" for 
the response function and the values of c^o scan the region of interest. In our case of triton 
photodisintegration we solve Eq. (0]) with F = 10 and 20 MeV and for a few hundred of ujq 
values chosen in the interval —40 MeV to 200 MeV. 
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Step 2. The overlaps of the solutions from Step 1. are calculated. Of course 

these overlaps depend on ujq and F. A theorem for integral transforms based on the closure 
property of the Hamiltonian eigenstates [? ] ensures that in the present case this dependence 
can be expressed as [? ] 

L(cJo,r) = (§1^) = J R{uj)C{uj,cuo,T)du;, (5) 
where C is the Lorentzian function centered at Uq and with F as a width: 

C{u,iJo,T) = - ^ . (6) 

- UJoY + ^ 

Therefore by solving Eq. one can easily obtain the Lorentz integral transform of the 
response function. 

Step 3. The transform (0) is inverted (see e.g. Ref. [? ? ]) in order to obtain the 
response function and with it the cross section. Of course the inversion result has to be 
independent of F and show a high degree of stability. 



III. THE HH EXPANSION TECHNIQUE 

In the HH expansion approach the basis functions for the solution of the A-body 
Schrodinger equation for the ground state |0), as well as of the A-hodj "LIT equation" 
for |\&), are expressed as functions of the hyperspherical coordinates. These are constructed 
starting from the A — 1 Jacobi coordinates C,i for the A-particle system. They consist of a hy- 



perradius pA-i = sJYlt ^(^)^ and a hyperangle VLa-i = {ii, 6, ■ ■ ■ , "^2, V^s, ■ ■ ■ , V^a-i} 
= ^A-i} , where the A — 2 hyperangles {{p2,(p3, ■ ■ ■ ,ipA-i} are defined in the 

[0,7r/2] interval through the relation 

sinv9„ = , n = 2,3, ■ ■ ■ , A - 1 . (7) 

In terms of these coordinates the Laplace operator for n Jacobi coordinates is 

A _ 1 ^ .3n-l d 1 ,;.2 



i.e. a sum of a hyperradial operator and of a term containing the hyperspherical, or grand 
angular momentum operator AT^. The latter can be expressed recursively in terms of the 
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hyperangular momentum operator associated with the subset {^i, • • • , Cn}) ^n-n ^n' 
the orbital momentum operator corresponding to the nth Jacobi coordinate, namely [? ]: 

^2 ^ ^ ^ 3n-6- (3n-2)cos293„ d 1 ^ 1 p rg^ 

" (9v92 sin2v9„ difncos^fn sinVn 

In the above equation Kf = i^, and the angular momentum operator associated with these 

— * — * — * 

n coordinates is L„ = L„_i + The operators K^, 1^, k^_i, I^i and V'n^ commute with 
each other. 

Given the set of hyperspherical coordinates, the hyperspherical harmonics functions y[Kn] 
are the eigenfunctions of the hyperangular momentum operator i^^, 

Kly[K^] = Kn{Kn + 3n - 2)y[K^] , (10) 

where Kn is the hyperangular quantum number and the symbol [Kn] stands for the set of 
quantum numbers {K^, X^-i, • • • , K2; L^, M„; -L„_i, • • • , -L2; • • • , £1}. The explicit 

expression for the HH functions of n Jacobi coordinates can be constructed recursivly as 
described in Ref. [? ]. 

The HH functions y[Kn] form a complete and orthonormal set of functions that satisfy 

{yiK^] I yiK^]) = SiK„],iK^] ■ (11) 

The particle permutation properties of the HH functions can be implemented in different 
ways. In the CHH approach used here this is done as described in Ref. [? ], while in the 
EIHH framework two recently developed powerful algorithms are employed [? ? ]• 

As functions of the hyperspherical coordinates the antisymmetric spin-isospin configura- 
tion space basis functions are written in the form 

^^KU7](P' ^A-i, a, t) = RNip)T^^l^^in, a, r) , (12) 

where, for the hyperradial part of the wave function one chooses an orthonormalized set of 
functions, Rn{p) ~ -^iv(p/Po) exp(— p/2po)- Here LJ^ arc the generalized Laguerre polyno- 
mials, and v and po are free parameters [? ]. In the fully antisymmetric functions r|}^^j 
the notation [K, 7] stands for the set of quantum numbers {K, J, Mj, L, 5", T, Mr}, i.e. total 
angular momentum and its third component, orbital angular momentum, spin, total isospin 
and isospin projection of the A-body system, respectively and o", r indicate collectively the 
spin and isospin variables of the A nucleons. These functions are obtained by the coupling of 
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the HH functions y^J!^^ (Vl) with specified L and M of the representation [/] and the spin- 
isospin functions ©sl/gTA/r ''") ^i^h specified S,Ms,T and Mt of the conjugate Young 
scheme [/]. In many cases the rate of convergence of the HH expansion is rather slow. In 
order to accelerate it one can follow two rather different approaches, leading to the CHH 
and EIHH methods. 



A. The CHH Approach 

The CHH approach [? ] has been used extensively in the past. While the underlying 
idea is common, there exist a number of slightly different versions of this approach. In the 
following we describe how the CHH has been applied by the present authors in this and 
previous works [? ? ]• 

The main idea of the CHH approach consists of the insertion in Eq. (fT^ of a Jastrow 
factor J' embodying the short range correlation due to the repulsive part of the potential. 
Such a repulsion leads to high momentum components in the wave function which make the 
convergence of the HH expansion rather slow. Therefore one modifies Eq. (|T^ into 

Cnn^ff^^Jp,nA-i,a,T) = J RNip)Tfl^p,a,T) . (13) 

Here JT" is taken as a state-dependent correlation operator of the form 

i<j s,t 

where Ps^t{ij) are projection operators onto nucleon pairs {ij) with spin s and isospin t and 
iS is a particle symmetrization operator. The radial parts of the correlation operator are 
chosen as follows. For rij < tq (the healing distance) fs,t{fij) is chosen to be the zero energy 
pair wave function in the corresponding s, t state. The healing is insured by imposing the 
condition fs,t{.fij) = 1 for r^j > Tq and j(rjj) = for r^j = r^. The s,t = 1,3 and s,t = 3,1 
cases are determined from the ^5*0 and ^Si partial wave of the NN potential. Since for the 
s,t = 1,1 and s,t = 3,3 cases the ^Pi and ^Pi potentials are not sufficiently attractive to 
obtain a healing distance, we introduce an additional intermediate range central interaction. 
Further details can be found in [? ]. 
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B. The EIHH approach 



In contrast to the CHH approach, which focuses on the short range properties of the wave 
functions, the EIHH approach focuses on the potential operator. The acceleration of the 
convergence is obtained by introducing a two-body effective interaction V^'^^ [? ? ] • This is 
done by the division of the total HH space in P and Q spaces realized via the HH quantum 
number K (P space: K < Kmax, Q space: K > K^ax )■ Since in hyperspherical coordinates 
the total Hamiltonian is written as 



H = i^[-\ + -r]+yv.,. (15) 




one can extract from it a "two-body" Hamiltonian of particles A and {A — 1) 



^2(P) = ^^ + I^A,(A-1), (16) 



which, however, contains the hyperspherical part of the A-body kinetic energy. Since the HH 
functions of the {A — 2) system are eigenfunctions of K\_2 one has an explicit dependence 
of H2 on the quantum number Ka-2 of the residual system, i.e. H2 —>■ H2^'^ . Applying 

the hermitian version of the Lee-Suzuki method [? ? ? ] to one gets an effective 

[21 

Hamiltonian -^26//- The effective interaction V^^^ is obtained from 

ySf'-\p) = Htff\P)-^^- (17) 



[91 _ _ 

This V^jy replaces Vij in Eq. ()15p when one projects the solution on the P-space. This effec- 

[21 

tive potential has the following property: V^^j — > Vij for P ^ 1. Due to the "effectiveness" 
of the operator the solution of the Schrodinger equation converges faster to the true one. 

The present HH formulation of an effective potential has the following peculiarities, which 

[21 

make it particularly "effective": (i) V^^j is p dependent, therefore it contains information 

on the "medium" ; (ii) because of the above mentioned Ka-2 dependence the (A-2) residual 

[2] 

system is not a pure spectator and its state influences the state of the pair, therefore V^j^ 
can be viewed as a state dependent effective interaction. Moreover we point out that the 
effective potential can be interpreted as a kind of momentum expansion, since the short 
range resolution is increased with growing K^ax- 
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IV. RESULTS 



We have compared the transform L(u;o,r) and response function R^u) as computed by 
the CHH and EIHH approaches. The potentials employed were the AV18 [? ] two-body 
and UIX [? ] three-body potential models. In the case of CHH we expand the ground state 
with 11 hyperradial functions and the hyperspherical states with K up to K^ax = 10, but 
consider additional selected states up to K^ax = 100. For the T=l/2 component of |\1') we 
take 26 hyperradial functions, K^ax = 9 and additional selected states up to Kmax = 65. 
For T=3/2 K^ax = 15 and 26 hyperradial functions are sufficient for a convergent result. 
In the case of EIHH we expand with 26 hyperradial functions up to Kmax = 14 (ground 
state) and K^ax = 15 (|^))- In Fig. [T]the EIHH rate of convergence of L(u;o,r) is shown 
for the (J=l/2,T=l/2) contribution to \I' with F = 10 MeV. One sees that the convergence 
is very good in the whole energy range. One also notes that small oscillations are present 
in the transform at higher energies. If needed one could improve the result by increasing 
the number of hyperradial states. Convergence patterns of the other (J,T) contributions 
in the EIHH calculation are similar to that shown in Fig. ^ Also the rate of convergence 
of the CHH calculation for T=3/2 is rather similar to that of Fig. ^ while for T=l/2 the 
convergence is considerably slower. Comparable differences in the rate of convergence of 
the two methods have also been found for the ^He total photoabsorption cross section with 
semirealistic potential models [? ]. Thus, in view of applications to more complex nuclei, 
the EIHH method seems to be preferable from the calculational point of view due to its 
faster convergence. In Fig. we show the results for the total Lorentz integral transform 
of R{u!) (contributions of J=l/2 and J=3/2 components are summed) for both the T=l/2 
and T=3/2 cases and with F = 20 MeV. It is evident that for both isospin channels the 
agreement is excellent over the whole considered energy range. In Fig. the relative 
difference AL = (L^™^ — L'-^^^)/L'-^^^ is shown in percent. It is readily seen that AL 
does not exceed 0.4% in the peak region. For higher energies the differences can become 
somewhat larger but remain below 2% up to pion threshold. In Fig. IHK we show the total 
photoabsorption cross section. The results of the two approaches remain very close also 
after inverting the transforms {Step 3 in Section |H}. In fact the EIHH and CHH curves 
are hardly distinguishable. Even if the agreement has slightly deteriorated compared to 
Fig. |2ldue to the inversion procedure, the quality remains excellent. As shown in Fig. Ob 
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the difference remains below 0.8% in the peak region and has a maximum of 2.4% up to 
pion threshold. We would hke to point out that the agreement of the two HH calculations 
in the peak region is even better than the already good agreement of CHH and Faddeev 
calculations shown in the above mentioned benchmark [? ]. We do not compare the present 
results with the Faddeev cross section in [? ] since that benchmark was made neglecting the 
charge dependent terms in the AVI 8 nuclear force. Though that effect is small (of the order 
of 1% in the peak region), it is much larger than the difference between CHH and EIHH 
calculations. If one considers that EIHH and CHH approaches differ substantially both from 
the conceptual and calculational points of view (of course completely independent codes are 
used), the agreement of the results is really remarkable and shows that one can rely on the 
application of the EIHH approach to heavier systems with a high degree of confidence. In 
particular, its application to the calculation of the total photonuclear cross section of ^He 
within realistic nuclar force models, containing two- and three-body forces, would be very 
much needed. It could shed some light on the role of three-nucleon forces, which might be 
even more important than already found in the three-nucleon system. Moreover, it would 
make a valuable contribution to the clarification of the present rather confused situation in 
the comparison between theoretical and experimental results in "^He (see e.g. [? ? ]). Work 
in this direction is in progress. 
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Fig.l 
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FIG. 1: Rate of EIHH convergence of the LIT of the unretarded dipole response function (Eq.(|SJ)) 
for r = 10 MeV, hmited to the (J=l/2,T=l/2) component. The even values of K^ax refer to the 
calculation of the ground state |0), the odd values to that of the Lorentz state ^. 
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FIG. 2: Comparison between CHH and EIHH results for the LIT of the unretarded dipole response 
function (Eq.©) in the two isospin channels, (a): CHH (full curve); EIHH (dashed curve), (b): 
The difference AL = (jj-^iiii -i^cnuy^cun between the CHH and EIHH results of Fig. 2a expressed 
in percent. 
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Fig. 3 
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FIG. 3: Comparison between CHH and EIHH total photonuclear cross sections of triton. (a): CHH 
(full curve); EIHH (dashed curve), (b): The difference AR = (ii!™" - i?C™)/^CHH between the 
CHH and EIHH results of Fig. 3a expressed in percent. 
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